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The abBA collaboration proposes to 
conduct a comprehensive program of 
precise measurements of neutron j3-decay 
coefficients a (the correlation between the 
neutrino momentum and the decay 
electron momentum), b (the electron 
energy spectral distortion term), A (the 
correlation between the neutron spin and 
the decay electron momentum), and B 
(the correlation between the neutron spin 
and the decay neutrino momentum) at a 
cold neutron beam facility. We have used 
a GEANT4-based code to simulate the 
propagation of decay electrons and protons 
in the electromagnetic spectrometer and 
study the energy and timing response of a 



pair of Silicon detectors. We used these 
results to examine systematic effects and 
find the uncertainties with which the 
physics parameters a, b, A, and B can 
be extracted from an over-determined 
experimental data set. 
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1. Introduction 

The abBA collaboration is proposing to perform a 
measurement of a "complete set" of correlations in the 
neutron j3-decay using the same apparatus, and improve 
the precision of the correlation coefficients a, b, A, and 
B by up to an order of magnitude. 

GEANT4 is a general-purpose software package for 
simulation of the passage of particles through matter 
that provides a complete set of tools for all domains of 
detector simulation [1]. In particular, the GEANT4 
toolkit currently provides particle tracking in non- 
uniform magnetic and electric fields and handles 
combined electromagnetic fields transparently [2]. The 
GEANT4 Low Energy Electromagnetic Physics group 
validates the low energy electromagnetic processes for 
electrons down to 250 eV [3]. 

In this report we describe a GEANT4 simulation of 
the abBA spectrometer and outline the algorithm for the 
extraction of the physics decay parameters. 



2. abBA Detector Geometry 

In the tentative design of the abBA spectrometer the 
decay particles (electrons and protons) are guided by 
the electric and magnetic fields and interact only in 
sensitive detectors thus avoiding the energy losses and 
scatterings in apertures, grids, or windows [4] . 

The simplified geometry of the detector defines 
two sensitive planar Silicon detectors with a 1 00 mm x 
100 mm^ area and a 2 mm thickness. The two Si detec- 
tors are separated by 4 m. The coordinate system is 
defined with the Si detectors at Xi2 = ± 2 m (Fig. 1). 

A passive solenoid magnet is placed around the 
decay region, with its axis of symmetry perpendicular 
to the incident neutron beam, along the x coordinate. 
The 3 m long magnet with a 0.8 m radius can produce 
a 4 T central magnetic field that decreases to 1 T at the 
detector positions, thus guiding charged particles from 
the decay region to the Si detectors. A tubular electrode 
held at ~ 30 kV accelerates the protons so they can be 
detected in a silicon detector. 
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Fig. 1. The schematic drawing of the main elements of the abBA spectrometer. 
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3. Magnetic Field 

The magnetic field along the x axis of the detector 
solenoid is given by: 

Y r ^ 



B(x,p=0) = 



2nNI 



R^+iL-xf 



VF77 



where L is the length of the solenoid, R its radius, 
and X the axial coordinate. Meanwhile, A'^ denotes the 
number of turns per unit length, and / is the electrical 
current in the closely wound cylindrical coil. 

Thanks to axial symmetry, the magnetic field off- 
axis, outside its sources, can be represented in terms of 
the magnetic field strength B(x, 0) along the axis: 

\2« 



3^ 
dx 



■■B{x)-^B\x,0)+- 



and 



'' '^' dp tl(n-iyn\ \2 



= ^5'(x,0)+* 



^+z^ 



where <^is the magnetic scalar potential and p = ^Jy 

is the axial radius coordinate. These fields have been 

programmed into the GEANT4 user routine. 

4. Event Generator 

The electrons from the neutron j3-decay are generat- 
ed from (5x5x5) mm^ central volume with the rela- 
tivistic differential decay rate given by [5]: 



dr 



dE^dn^dn^^ 



(O.V.,f F{E^)\P.\E. 



_ V^F'^ud 



{2ny mSE+E^+ESp-pS 



■\M\ 



where E^ and p^ (E^ and p^) are the electron (neutrino) 
energy and momentum, m^ is the neutron mass, Gp is the 
Fermi constant, V^^ is the Cabbibo-Kobayashi-Maskawa 
matrix element, P=pJE^, and F{E^ is the Fermi func- 
tion that describes the interaction of the electron and the 
recoil proton. 

The transition matrix element squared \Mf' is given by 



lA/r =«.„«.,£,£.! l+|-e^l+|-5« 



xC„(£J(l+3gl) 



{l+[l+|-.f 



c;iK)Pk 
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where m^ is the proton mass, a is the fine structure 
constant, e^ is a low energy constant, 5„'s are model- 
independent radiative corrections, gp^ is the axial coup- 
ling constant, and the correlation coefficients a. A, B 
are incorporated into the recoil corrections C/^'e) [^l- 

5. Results and Conclusions 

A GEANT4 simulation of abBA detector energy 
and timing response was performed for 10^ neutron 



j8-decays. We used the values of the correlation coeffi- 
cients from Ref [7] (a=-0.1039, Z? = 0, ^ = -0.1161, 
B = 0.9878). For each event we recorded the neutron 
polarization, generated momenta of the final state 
particles and measured energy depositions and timing 
hits in the Silicon detectors. 

The separate GEANT4 run which included systemat- 
ic effects (energy and timing resolutions of Si detectors, 
detector calibration uncertainties, detector response 
nonlinearities, magnetic field inhomogeneities, neutron 
polarization uncertainty, etc.) was used to simulate the 
experimental data. ("MC data"). The flow chart of the 
physics analysis is summarized in Fig. 2. The correla- 
tion coefficients and their fitted uncertainties are 
extracted using the standard MINUIT code [8]. 
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Fig, 2, A flow diagram of the analysis. Monte Carlo histograms are recalculated in each 
step of MINUIT minimization. The "experiment" depends only on measured quantities 
(subscript £^, while the MC column depends both on generated variables (subscript 7) and 
simulated detector response (subscript £). 
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We present two example results: (i) the extraction of 
the parameter b with unpolarized neutron beam in 
Fig. 3, and (ii) the asymmetry coefficient A for 80 % 
polarized neutron beam in Fig. 4. At the current stage 
of development, a GEANT4 simulation limited to 1 0^ 
neutron decay events and 10^ MC data events, runs 
24 CPU hours on a 1 GHz Linux computer. Given the 



limited event statistics, analysis of MC data results in 
the coefficient b = -0.0025 ± 0.0028 and the coefficient 
A =-0.1170 ± 0.0010, where statistical and systematic 
uncertainties are combined. The code will be made 
faster by using the adiabatic invariants for charged 
particle tracking in the electromagnetic field, which 
will markedly improve the uncertainties of our method. 
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Fig. 3. Monte Carlo energy spectrum of neutron ^decay electron (top) and fractional differences 
A = (exp - the)/(exp + the) between the simulated and "MC data" spectrum (bottom). 
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Fig. 4. Monte Carlo histogram of the p asymmetry (top) and fractional differences between the 
simulated and "MC data" spectrum as a function of the electron energy (bottom). 
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